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Parametrically-excited surface waves, forced by a repeating sequence of -/V delta-function impulses, 
are considered within the framework of the Zhang- Vihals model 1|. The exact impulsive-forcing 
results, in the linear and weakly nonlinear regimes, are compared with results obtained numerically 
for sinusoidal and multifrequency forcing. We find surprisingly good agreement between impulsive 
forcing results and those obtained using a two-term truncated Fourier series representation of the 
impulsive forcing function. With impulsive forcing, the linear stability analysis can be carried out 
exactly and leads to an implicit equation for the neutral stability curve. As noted previously by 
Bechhoefer and Johnson 0, in the simplest case of N = 2 equally-spaced impulses per period (which 
alternate up and down) there are only subharmonic modes of instability. The familiar situation 
of alternating subharmonic and harmonic resonance tongues emerges only if an asymmetry in the 
spacing between the impulses is introduced. We extend the linear analysis for N = 2 impulses 
per period to the weakly nonlinear regime, where we determine the leading order nonlinear satura- 
tion of one- dimensional standing waves as a function of forcing strength. Specifically, an analytic 
expression for the cubic Landau coefficient in the bifurcation equation is derived as a function of 
the dimensionless spacing between the two impulses and the fluid parameters that appear in the 
Zhang- Vinals model. As the capillary parameter is varied, one finds a parameter region of wave 
amplitude suppression, which is due to a familiar 1:2 spatio-temporal resonance between the sub- 
harmonic mode of instability and a damped harmonic mode. This resonance occurs for impulsive 
forcing even when harmonic resonance tongues are absent from the neutral stability curve. The 
strength of this resonance feature can be tuned by varying the spacing between the impulses. This 
finding is interpreted in terms of a recent symmetry-based analysis of multifrequency forced Faraday 
waves lij-LJ 

PACS numbers: 05.45.-a, 47.35. +i, 47.54.+r, 89.75.Kd 



I. I. INTRODUCTION 

Standing waves form spontaneously on the free surface 
of a fluid layer when it is subjected to a time-periodic ver- 
tical vibration of sufficient strength. The waves, which 
result from a symmetry-breaking parametric instability, 
organize themselves into remarkably regular patterns, as 
first described by Faraday in 1831 In modern ex- 
periments, a variety of container geometries, fluids and 
forcing functions have been employed, resulting in a rich 
variety of standing wave patterns, as well as a range of 
dynamical responses. See 0, Q for reviews on Faraday 
waves. In this paper, we investigate the response of the 
Faraday system to an idealized forcing function consist- 
ing of a periodic sequence of delta- function impulses. The 
results are contrasted with those obtained in the classical 
cases of single and multifrequency forcing. 

The versatility of the Faraday system for investigat- 
ing pattern formation owes much to the vastness of its 
control parameter space. By varying the forcing fre- 
quency and amplitude, as well as the fluid properties, ex- 
perimentalists employing a sinusoidal acceleration func- 
tion have teased out spatially regular patterns ranging 
from stripes, hexagons and squares to targets, spirals 
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and quasipatterns d, 0, 0, ITU . The seminal experi- 
ments of Edwards and Fauve [12| showed how the addi- 
tion of a second commensurate frequency component in 
the forcing function could lead to even greater versatility 
of this system, via the controlled introduction of addi- 
tional length scales. Two-frequency forcing has led to a 
variety of exotic superlattice patterns 0, fidl . as well as 
quasipatterns Il5lll6l , triangular patterns |l7j and local- 
ized structures [l8l |. A detailed description of patterns 
readily achieved in two-frequency Faraday experiments 
can be found in |19| . Recent experiments by Epstein and 
Fineberg poj have shown how a third perturbing fre- 
quency can be used to rapidly switch between the novel 
patterns achieved with two-frequency forcing. 

Theoretically, the Faraday system presents a number 
of challenges due both to the explicit time-dependence of 
the forcing and the free-boundary nature of the problem; 
for a partial review of theoretical work see [3- Con- 
sequently, models incorporating simplifying assumptions 
have often been relied on, in conjunction with numeri- 
cal linear stability analysis or perturbative methods. In 
this manner, much has been accomplished. Linear sta- 
bility results build on the classic paper of Benjamin and 
Ursell [2l| , who explained Faraday's observation of a sub- 
harmonic response of the fluid layer to sinusoidal forcing. 
Specifically, they showed that the linear stability prob- 
lem, in the case of an ideal fluid, reduces to a Mathieu 
equation, for which the natural frequency of oscillation 
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is given by the inviscid dispersion relation. Subsequent 
investigators carried out linear stability analyses on the 
full fluid problem with sinusoidal forcing HI S E3 • 
These results were also extended to two- frequency forcing 
in pr| . Theoretical understanding of the nonlinear prob- 
lem progressed with the introduction of a quasipoten- 
tial formulation of the Faraday wave problem by Zhang 
and Vinals |l|. This model, which describes small am- 
plitude surface waves on a semi-infinite, weakly viscous 
fluid layer, eliminates the bulk flow, assumed to be poten- 
tial, from the description. The free boundary z = ft,(x, t) 
(x ? I 2 ) is then prescribed by an evolution equation for 
h, which is coupled to an evolution equation for a surface 
velocity potential see Section |n] Zhang and Vinals 
used this model to further investigate nonlinear effects 
via an asymptotic expansion, which demonstrated the 
importance of resonant triads (three-wave resonance ) in 
the pattern selection process. Later, Chen and Vinals |27j 
used perturbation methods to derive amplitude equations 
directly from the Navier-Stokes equations assuming only 
that the depth of the fluid layer is (effectively) infinite. 
These approaches yield good agreement with experiments 
in the relevant parameter regimes, but typically require 
a numerical evaluation of relevant quantities to make the 
comparison. 

An alternative theoretical approach to the Faraday 
wave pattern selection problem is based in equivariant 
bifurcation theory |28l |29| . This "model- independent" 
treatment of the symmetry-breaking bifurcations rele- 
vant to the Faraday experiment has proved especially 
fruitful in analyzing spatially-periodic superlattice pat- 
terns |33, 0S1 • An extension of this approach, tailored 
to the multifrequency Faraday problem, uses spatio- 
temporal symmetry arguments to derive amplitude equa- 
tions based on the patterns' symmetries. The form of the 
coefficients in the equations reflect the temporal (param- 
eter) symmetries of the forcing function and the time- 
reversal symmetry that is broken by the damping. The 
equations describe not only the critical modes, but also 
certain damped modes that are nonlinearly driven by 
them. These studies reveal the role of the various forcing 
parameters (frequencies, amplitudes, and relative phases) 
in nonlinearly coupling the driven modes to key weakly- 
damped modes. The results give insight into the pat- 
tern selection process and can be used to explain a num- 
ber of characteristics of the experimentally observed pat- 
terns Hum imhi. 

The investigations reviewed above focus on Faraday 
waves parametrically-excited by sinusoidal and multifre- 
quency forcing functions. An alternative forcing func- 
tion, consisting of a periodic sequence of delta-function 
impulses, was proposed by Bechhoefer and Johnson Q 
who indicated how the linear analysis could be greatly 
simplified in this case. (This idealization of parametri- 
cally forced systems has been put forth in a variety of 
physical contexts; see, for example, [3^, |3|| 137} .) Re- 
cently, Huepe and Silber j3g| showed that the linear anal- 
ysis of the impulsively-forced Faraday problem for a vis- 



cous fluid, as proposed in J2j, breaks down in certain 
parameter regimes. In these regimes, the flow in the 
fluid bulk cannot be (instantaneously) matched across 
the delta impulses. However, these complications are ab- 
sent from the Zhang- Vinals formulation, which describes 
the evolution of the free surface, with the bulk flow elim- 
inated from the model. In this paper we extend the re- 
sults of Bechhoefer and Johnson on impulsively-forced 
Faraday waves into the weakly nonlinear regime, within 
the framework of the Zhang- Vinals model. 

Following the method of |3^| , we develop a simple mod- 
ular approach to constructing the stroboscopic map asso- 
ciated with the linear stability problem. Specifically, we 
construct the linear maps which evolve perturbations of 
wavenumber k from one impulse in the sequence to the 
next. The linear map associated with one period of the 
forcing is then obtained by multiplication of N matrices, 
one for each impulse in the periodic sequence. From this 
simple construction we can derive explicit expressions for 
the Floquet multipliers associated with the linear stabil- 
ity problem, and arrive at an equation (implicit if N > 2) 
that describes the neutral stability curve. We analyse the 
results in the simplest case of N = 2 impulses in detail, 
and compare them with those obtained with an M-term 
truncated Fourier series representation of the impulsive 
forcing function. While the neutral stability curves are 
poorly approximated by the truncated Fourier series even 
for very large M, the onset parameters of critical forc- 
ing strength and wavenumber are quite accurately deter- 
mined with a single-term sinusoidal approximation (i.e., 
for M = 1). 

We extend the analysis for N = 2 impulses per period 
to the weakly nonlinear regime for one-dimensional sur- 
face waves. A key observation of Bechhoefer and John- 
son 0] is that there are no harmonic resonance tongues 
for a sequence of N = 2 impulses if they are equally 
spaced in time (i.e., an up impulse and then a down im- 
pulse half a period later). We show that, despite the 
absence of harmonic instabilities, there can still be a res- 
onant interaction between the excited subharmonic and 
damped harmonic modes which leads, in one dimension, 
to an associated degradation in the nonlinear response of 
the fluid to the vibration in the weakly nonlinear regime. 
This is similar to the spatio-temporal 1:2 resonance that 
occurs in the case of sinusoidal forcing Q] , for which sub- 
harmonic and harmonic instabilities alternate with in- 
creasing perturbation wavenumber k. Thus, although 
the harmonic modes can never be linearly excited in this 
particular impulsive case, these modes can nonetheless 
be driven nonlinearly by the critical subharmonic modes 
and in this way influence the nonlinear evolution of the 
instability. In addition to analysing the case of N = 2 
equally- spaced impulses, we explore the effects of vary- 
ing the spacing between the impulses so that they are 
no longer exactly half a period apart. We find that this 
asymmetry in impulse timing enhances or diminishes the 
influence of the 1:2 resonance on the standing wave am- 
plitude depending on the sense in which it is applied. We 
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understand this result by considering a two-term Fourier 
series approximation to the impulsive forcing function, 
and applying recent results pertaining to multifrequency 
forcing 0,U- 

Our paper is organized as follows. In the next section 
we present the Zhang- Vihals model of Faraday waves, 
and then carry out the linear stability analysis in the 
general case of N impulses per forcing period. We then 
focus on the simplest case of N = 2 and examine in 
detail the sequence of subharmonic and harmonic reso- 
nance tongues as a function of the dimensionless spac- 
ing between the impulses. We compare our results with 
those obtained with sinusoidal and multifrequency forc- 
ing functions within the context of a truncated Fourier 
series representation of the impulsive forcing function. In 
section III we derive the cubic bifurcation equation that 
determines the amplitude of one-dimensional spatially- 
periodic surface waves driven by impulsive forcing, fo- 
cusing on the 1:2 spatio-temporal resonance feature. We 
compare our weakly nonlinear results with those obtained 
with single and two-frequency forced Faraday waves. Fi- 
nally, in Section IV, we summarize our results and indi- 
cate some promising directions for subsequent investiga- 
tions. 



II. LINEAR RESULTS FOR THE 
ZHANG- VINALS FARADAY WAVE MODEL 



This section introduces the Zhang- Vihals model for 
Faraday waves. A linear stability analysis is then per- 
formed for waves excited by a periodic sequence of N 
impulses and detailed results are presented for the sim- 
plest case of iV = 2. 



A. Zhang- Vinals model 

The quasipotential formulation of the Faraday wave 
problem, due to Zhang and Vihals 0, is derived from 
the Navier-Stokes equations assuming small amplitude 
surface waves on a deep, nearly inviscid fluid layer. It de- 
scribes the free surface height h(x, t) and surface velocity 
potential <&(x, t) of a fluid subjected to a (dimensionless) 
periodic vertical acceleration function G(t). Employing 
these equations greatly simplifies our calculations since 
the flow in the bulk does not appear explicitly — we 
need only track the behavior of the free surface, a func- 
tion of the horizontal coordinate x € R 2 and time t. The 
Zhang- Vinals model takes the form 



where the nonlinear terms in Q are given by 
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Here the operator D multiplies each Fourier component 
by its wavenumber, e.g. De = |k|e l , and the di- 
mensionless fluid parameters are 
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where g is the usual gravitational acceleration, lo is the 
forcing frequency (27r divided by the forcing period), v 
is the kinematic viscosity, p is the density, and T is the 
surface tension. The wavenumber k$ is chosen to satisfy 
the inviscid dispersion relation 
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p 



(4) 



where % is the frequency associated with the typical sub- 
harmonic response of Faraday waves. After dividing (J3J 
by u> 2 we find Go + To = 4. G(t) describes the applied 
acceleration; since time has been scaled by the forcing 
frequency ui, the dimensionless period of G(t) is 2it. 
We write the impulsive forcing function in the form 



Gi mp (t) — fn${t ~ tfi) 



(5) 



It is parameterized by the locations t n of the impulses 
{t n < tn+i, t n +N — t n = 27r) and by the amplitudes f n 
(fn+N = fn), where N is the number of impulses in the 
repeating sequence. The dimensionless amplitude /„ of 
an impulse is given by 

, _ v n k 

fn = , (6) 

where v n is the jump in velocity at time t n . We must 
also require 
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0, 



(7) 



to prevent a net translation of the container (for suitably 
chosen initial velocity). An example of an impulsive ac- 
celeration function, along with the corresponding velocity 
and position functions, is shown in Fig. ^ 



B. Linear stability analysis 

1. Calculations 

Our linear stability calculations follow the method 
of |2|. From (JIJ, the linear stability problem takes the 
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FIG. 1: Example of (a) acceleration, (b) velocity, and (c) 
position of the fluid container over time in the case of two 
impulses per period. The positions of the vertical lines in 
the acceleration function denote the locations of the delta 
functions; here they have equal magnitude. 



form 

\d t - 7 V 2 ) 2 -^(r V 2 -(Go-G(i)))] h{x,t) = 0. (8) 

We then write h(x,t) — pfe(t)e lkx + c.c., where c.c. de- 
notes complex conjugate, and find that Pk(t) satisfies 

p'l + 2 7 fcV fe + [7 2 fc 4 + Tofc 3 + (Go - G{t))k]p k = 0, (9) 

where the prime denotes differentiation with respect to 
t and the k subscript emphasizes the dependence of the 
solution on the perturbation wavenumber k — |k|. Here- 
after we focus on k ^ 0; the k — mode cannot be 
excited due to mass conservation. (This is reflected in 
the governing equations J2J) through the spatial deriva- 
tives in the nonlinear terms which ensure that the k = 
mode decouples from the other modes.) 

Between each impulse, (Q is simply the equation for a 
damped harmonic oscillator, with solution 



Pk 



(t) = A k , n e^ k2+l ^ t - t ^ +c.c, t€(t n ,t n+1 ). (10) 



Here uj k is the natural frequency of a wave with wavenum- 
ber k, which is determined by the dispersion relation 



r fc 3 + G k. 



(ii) 



We demand that Pk(t) be continuous across each im- 
pulse, i.e., p k (t+) = Pk{t~) = p k (t n ), where p k (tt) = 
lim t ^ t ± Pk(t). Integrating 10 across the n th impulse we 
obtain the following jump condition: 

p' k (4)-p' k (t-) = fnk Pk (t n ). (12) 

This condition, together with the continuity requirement, 
yields the following map from A ky7l to A kt . n+ i: 



k,n+l j =e -~(k*d nMi 
k,n+l 



k . n 



where 

M kn = 



Sk, 



Ck,r< 
F k ' 



-lCfc,: 



C ky n + F kyn+ iS k . 



(13) 
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Here A r k n (A], n ) is the real (imaginary) part of A kjTl 
and c kin = cos(u> k d n ), s fc ,n = sin(w fc d„), F ky7l = 
Note that l|14|l depends on only two parameters: d n = 
tn+i—tn: the time between the n th and (n+l) st impulses, 
and fn+i, the strength of the (n + l) s * impulse. 

Piecing together the relationships between A n and 
A n+ i, A n+ i and A n+ 2, etc. across one period, we find 
the stroboscopic map that relates the solution of the lin- 
earized problem at time t n to the solution one period 
later (i.e., after the sequence of N impulses): 
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M k)U , where M k _j is given 



Here M k 
by IH|). 

The eigenvalues of e~ 2 ^ lk M k determine the linear 
stability of the interface to disturbances of wavenum- 
ber k. Note that M k has determinant one, so these 
eigenvalues (A±) can be related to the trace of M k by 



e 27r 7 fc A± = iTr(M fc ) ± ^^Tv(M k )) 2 - 1. An instabil- 
ity exists when both eigenvalues (Floquet multipliers) are 
real and the magnitude of one exceeds unity. The thresh- 
old condition is therefore 
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Tr(M fc ) = ±cosh(27T7fc 2 ), 



(16) 



where '+' corresponds to the harmonic case (Floquet 
multiplier +1), and '— ' corresponds to the subharmonic 
case (Floquet multiplier —1). This threshold condition 
defines an implicit equation for the neutral stability curve 
/(fc), where / is an appropriate measure of the overall im- 
pulse strength (e.g., for TV = 2 we use / = |/„| = |/ n+ iJ). 

Examples of neutral stability curves for various acceler- 
ation functions are presented in Figs. l2l4l -l5l For values of 
/ below the minimum of these curves, the trivial solution 
(h = $ = 0) is stable. When / is increased above the crit- 
ical forcing amplitude f c at the curves' minimum (with 
corresponding critical wavenumber k c ), the flat surface 
solution becomes unstable to standing waves. An exam- 
ple of the associated critical eigenfunction, p k (t) (k = k c 
and / = f c ), is shown in Fig. |5J Note the presence of 
kinks in p k (t) at the impulses, a consequence of the jump 
condition (|12|l . 



2. Linear stability results for N — 2 

We now present a detailed analysis for the case of two 
impulses per period. In this illustrative example, the 
forcing function is parameterized by two constants: /, 
which characterizes the magnitude of the velocity jump 
associated with each impulse (see @)), and 



A 



dp - dx 
2tt 



e(-l,i) 



(17) 
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which measures the deviation from equal spacing of the 
impulses. Thus A = corresponds to equally-spaced im- 
pulses, and as |A| increases the spacing between the im- 
pulses becomes increasingly asymmetric; see Fig.^for an 
example where A < 0. Specifically, we set /„ = (—1)™/, 
and 



in, — 



7r(n + (l-(-l)*)^) (18) by 



in lp|l. where we have assumed, without loss of generality, 
that t = 0. 



It follows from 115fl that Ak :H +2 and Ak, n are related 



where 



A k,n 



A k,n+2 \ _ -27r 7 fc 2 ( COs(27TWfc) - F n C n S n+1 - sm(2nuj k ) + F n 

A k,n+2 ) ~ V sin(27Tw fc ) + F n s n+ i(s n + F n c n ) cos(27rw fc ) + F n s n+ i(c n - F„s n ) J \ A\ n 



(19) 



fnk 

c„ = cos(w fc (i„), s n = sin(u; k d n ), d n = n(l + (-1)™A), F n = . (20) 



Solving (|16|) for the neutral stability curve /(fc), we find 

f(]A - 1^1 / _ cosn ( 27r 7 fc2 ) ~ cos(27ra; fc ) 
^ ' ~ k y cos(27rw fc ) - cos(27rw fe A) ' ^ ' 

where ujk is given by (|llf> . As before, '+' corresponds to 
the harmonic case and '-' to the subharmonic case. Note 
that f(k) is invariant under A — > —A, so we can restrict 
to A e [0, 1) for the remainder of the linear analysis. This 
is not, however, a symmetry of the full problem and we 
will see that the sign of A affects the nonlinear response 
of the fluid to the impulsive forcing function. 

The expression for the neutral stability curve /(fc), 
given by (|2~T1) . diverges at each value of k for which 
cos(27raj fe ) = cos(27ru; fc A). This leads to the structure 
of clearly demarcated resonance tongues seen in Fig. [21 
The presence of vertical boundaries between resonance 
tongues contrasts sharply with the classical case of sinu- 
soidal forcing (Fig.0{a)), where multiple instabilities (in- 
volving distinct tongues) can occur for a single wavcnum- 
ber. We now examine in more detail the spacing and se- 
quence of resonant tongues in the case of two impulses per 
period. It proves useful, for the purposes of this analysis, 
to use LUk as the independent variable, where Uk and k are 
related through the monotonic dispersion relation I jllj l. 

The divergence of the neutral curves, which leads to 
vertical asymptotes separating the resonance tongues, is 
due to the denominator in expression l|21|l . We therefore 
introduce 

D(ujk) = cos(27ro;fc) — cos(2iTUJk A), (22) 

and denote the successive zeroes of D(ujk) by uj 3 k (u> 3 k < 
to J k +1 , uj k = 0). Since for k > the numerator of the 
expression inside the square root in H21J) is strictly pos- 
itive in the harmonic case and negative in the subhar- 
monic case, it follows that D(u>k) > for harmonic reso- 




FIG. 2: Examples of neutral stability curves f(k) from I21H 
for various two-impulse forcing functions. Solid curves corre- 
spond to subharmonic tongues (Floquet multiplier —1), while 
dashed curves indicate harmonic tongues (Floquet multiplier 
+ 1). The spacing of the impulses, given by I17L is (a) A = 0, 
(b) A = 0.05, (c) A = i, (d) A = |. The dimensionless 
parameters in the calculations are 7 = 0.02 and To = 0.04 
(corresponding to fluid parameters v = 0.1 cm 2 /s, r = 16 
dyn/cm, p = 1 g/cm 3 , and ^ = 20 Hz). 



nance tongues and D(ujk) < for subharmonic resonance 
tongues. 

In the general case, resonance tongues will alternate 
between harmonic and subharmonic as k (equivalently, 
u>k) increases and D(u>k) alternates in sign. However, 
for special values of A, successive subharmonic (or har- 
monic) tongues may occur. Specifically, this happens for 
some j if -D(^) = Z?'(w^,) = 0. We denote these double 
zeroes by c2>™, m € Z + (a)™ < w[ n+1 ). It is straight- 
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FIG. 3: Eigenfunction Pk(t) at k = k c , f = fc for an impulsive 
forcing function with equally-spaced impulses. The dimen- 
sionless parameters used are 7 = 0.17, To = 0.19 (roughly 
corresponding to physical parameters v = 0.20 cm 2 /s, Y = 16 
dyn/cm, p = 1 g/cm 3 , and — 80 H). Note the disconti- 
nuities in the derivative of Pk(t) induced by the impulses at 
t = nir. 



forward to show that these degenerate values occur only 
when cos(27rw™) = cos(27rd;™A) = ±1, which implies 
A = I for some co-prime integers p, q (0 < p < q) . 
There are thus two possibilities: either both p and q are 
odd, or one of them is even. If they are both odd, then 
= successive subharmonic (harmonic) tongues 
occur when m is even (odd). (See Fig. ^c), generated 
with A = -|.) In contrast, if either p or q is even (or 
p = 0), then only subharmonic tongues can occur in 
succession; these are separated by vertical asymptotes 
at cj™ = mq. (See Fig. |3[a) and (d) for examples of 
this case.) Note that when A = there are vertical 
asymptotes at each integer value of Wk separating succes- 
sive subharmonic tongues. In other words, there are no 
harmonic resonance tongues for equally-spaced impulses 
(Fig-GIa)), as previously noted by Bechhoefer and John- 
son 

The maximum possible width for the resonance 



i+i . 



,3 - 



1. 



tongues occurs when A = 0, in which case ui k 
As A is perturbed away from 0, the vertical asymptotes 



(located at Loi 



j for A = 0) split, and a thin 



harmonic tongue appears between them thereby separat- 
ing the two fat subharmonic tongues (see, for example, 
Fig-[2Ib)). Analogous results hold when A is perturbed 
away from other rational numbers. 



3. Comparison with sinusoidal and multifrequency forcing 

We now compare the linear stability results for equally- 
spaced impulsive forcing (i.e., A = 0) with the corre- 
sponding results for sinusoidal forcing 



Gsin(i) = /sinSin(t). 



(23) 



Here / S j n = g3 ^ fc ° , where <? S i n is the maximum (dimen- 
sioned) acceleration. (See Fig. 0Ja) for an example of a 
neutral curve in the case of sinusoidal forcing.) As noted 
earlier, the most striking difference between these two 
cases, manifest in the neutral stability curve ipTJl. is the 





FIG. 4: Neutral stability curves for truncated Fourier series 
approximations of equally spaced impulses with (a) 1 term, 
(b) 2 terms, (c) 10 terms, and (d) 20 terms. Large points 
correspond to subharmonic tongues and the smaller points 
to harmonic tongues. Parameters are the same as those in 
Fig. HI C urves were generated using the method described 
in gall- 




FIG. 5: Neutral stability curve for impulsive forcing with 
A = (solid lines) and a truncated Fourier series approxi- 
mation with 50 terms (dots). Note that the 5th (rightmost) 
resonance tongue is not captured at all by the approximation. 
Parameters are as in Fig. [5] 



lack of harmonic tongues for the impulsive case. This dif- 
ference is a consequence of how the Floquet multipliers 
move in the complex plane as the perturbation wavcnum- 
ber k is varied, as illustrated in Fig. In both the si- 
nusoidal and the impulsive cases, the Floquet multipliers 
at fc = (labeled as point 1 in Figs. Efa,b)) are degen- 
erate and take the value +1. As fc is increased from 0, 
these split into a complex conjugate pair lying inside the 
unit circle which, together, sweep out a path encircling 
the origin. When fc reaches a certain value the Floquet 
multipliers meet again on the negative real axis and, at 
this point (labeled 2), are degenerate. With additional 
increase in fc, they separate and move in opposite direc- 
tions along the (negative) real axis. If the magnitude of / 
exceeds the critical value (as it does in Figs.[IJa,b)), then 
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FIG. 6: Floquet multipliers fi± of Eq. 1241 plotted in the com- 
plex plane over the range k £ (0, 3] with / = f c (l + 1.4) (this 
forcing is well above the onset value for the first resonance 
tongue but below onset for the second) for (a) equally-spaced 
impulses, and (b) sinusoidal forcing. These cases differ in that 
only for sinusoidal forcing do the Floquet multipliers split on 
the positive real axis. Fluid parameters are as in Fig.0 Num- 
bers denote points described in the text. 
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FIG. 7: (a) Ratio of critical forcing strength for impulsive 
and sinusoidal forcing as a function of forcing frequency. Vis- 
cosities are labeled in the figure; other parameters are as in 
Fig. [5] The solid line for v = lcS is indistinguishable from 
a horizontal line at ?. (b) Log-log plot of 



11 _ lU 



T fc,sin 1 

function of 7 with To = 0.13 (solid line). For comparison, a 
line of slope two is drawn nearby (dashed line). 



the most negative Floquet multiplier eventually crosses 
— 1, indicating a subharmonic instability. In any case, 
as k is increased further, these (real) Floquet multipli- 
ers soon reverse direction and recombine (point 3) before 
splitting into a complex conjugate pair that moves to the 
right, again encircling the origin, before merging on the 
positive real axis (point 4). It is here that the difference 
between the two cases emerges. With sinusoidal forc- 
ing the Floquet multipliers split and move along the real 
axis; a harmonic instability will therefore develop if the 
forcing is strong enough to drive the rightmost Floquet 
multiplier outside the unit circle. In contrast, the Flo- 
quet multipliers in the impulsive case do not split along 
the positive real axis but simply continue their trajectory 
around the origin as a complex conjugate pair. 

The above description follows directly from an exami- 
nation of the Floquet multipliers in (|19fl : 



M± = e - 2 ^ fc2 (a ± \/«2 _ 1) ( (24) 



where, for A = 0, 

1 ^ / , , . if fk sin(7raji,) \ 2 

a = -Tr M fe = cos(2™ fe ) - - ± < 1. 

2 2 \ ojk J 

(25) 

From l|24(l we see that the Floquet multipliers are real 
and distinct whenever a 2 > 1. However, since a < 1, 
this can happen only when a < —1, in which case they 
are negative. 

To further compare equally-spaced impulses and sinu- 
soidal forcing, we consider the respective critical forc- 
ing strengths at the onset of standing waves, f c ,imp a- n d 
fc,sin- In the limit of weak damping we may expand l|21|) 
about 7 = and ui^ = ^ (near the minimum of the neu- 
tral stability curves for small 7) to obtain 

fc.zmp = 7rfc c7 (l + 0(7 2 ))- (26) 
For sinusoidal forcing and small damping we have 

/ c , s *„ = 2fc c7 (l + 0(7 2 )), (27) 

(see, for example, 0). Thus f -ff^ = § + 0(7 2 ). This 
is the same ratio that arises when considering (the first 
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FIG. 8: Ratio of critical forcing strengths for equally-spaced 
impulses and the M-term truncated Fourier series, as a func- 
tion of M, for T = 0.04 and 7 = 0.02 (filled circles), 7 = 0.04 
(triangles), 7 = 0.1 (open circles). 



term in) the Fourier series expansion of the two-term 
impulsive forcing function with A = 0: 



G tmp (t) = -f imp V cos((2j + l)t). (28) 

7T * — ' 

j=0 



Here we have used (0 and (|18|l with /„ = (— l) n fi m p- ln 
Fig. EJa) we plot the ratio ^ C| "" p as a function of forcing 

frequency (the most easily tuned parameter in experi- 
ments) for several viscosities. As anticipated, for small v 
this ratio is well approximated by \. Furthermore, the 
ratio decreases with increasing viscosity (equivalently, 7) 

with a deviation |1 - ^ff^\ = 0{~f 2 ) (see Fig.Efb)). 

We may also consider the effect of adding a second 
term, cos(3i), to the truncation of the Fourier series l)28|) 
and examine / C) 2, the critical value of this two- frequency 
forcing function, in the weak damping limit. (In gen- 
eral, fc.M denotes the dimensionless forcing strength at 
onset for the M-term truncation of It is demon- 

strated in j3^| that for a two-frequency forcing function 
with commensurate frequencies mu> and ntu (muj is as- 
sumed to drive the primary instability), the nui frequency 
component is destabilizing when — < \/2 w 1.19. By this 
we mean that the threshold for instability is lower in the 
two-frequency case than in the single frequency case. For 
the two-term truncation of l|28[l , 21 = 1 and hence we ex- 
pect that / Cj2 < / c ,i = § fc,sin, i-e., jf^- < f • Indeed, 
this is consistent with our finding that y'"" p < ^ for 
impulsive forcing (see Fig.[7{a)). 

As expected, adding more terms to the truncation of 
the Fourier series l|28|) results in a closer approximation 
to the critical forcing strength f c ,imp, as shown in Fig. IH1 
Despite this good agreement between f c ,i m p an d f c ,M, 
the M-term truncated Fourier series is far less successful 
in capturing the behavior of the resonance tongues away 
from the minimum of the first tongue; this can be seen 
in Fig. El For 7 = 0.02 and r = 0.04, for example, / C)Af 



approximates f c ,imp remarkably well even for small M 
(filled circles in Fig.lHl; however, with M = 50, only the 
first four resonance tongues plausibly resemble those for 
impulsive forcing (see Fig. |SJ. 



III. WEAKLY NONLINEAR ANALYSIS 

A. Weakly nonlinear calculation 

We now extend our analysis of N — 2 impulses to the 
weakly nonlinear regime in the case of one-dimensional 
waves. 

Following we perform a multi-scale expansion: 

h(x, t, T) = ehx(x, t, T) + e 2 h 2 (x, t, T) + ■ ■ ■ , (29a) 

$(x,t,T) = e$ 1 (a;,t,T) + e 2 $ 2 (iM,T) + -- - , (29b) 

where e -C 1, the slow time is T = e 2 t, and the forcing 
amplitude is / = / c (l + e 2 /2)- We seek spatially-periodic 
solutions in separable Floquet-Fourier form: 



hi (x, t, T) = Z{T) Pl (t)e ik ° x +c.c, 
h 2 (x, t, T) = Z 2 (T) P2 (t)e 2lk ^ + c.c, 
$i(x,t,T) = Z(T) qi (t)e^ x + c.c, 
<S> 2 {x 7 t,T) = Z 2 {T)q 2 (t) 



,2ik r 



+ C.C, 



(30a) 
(30b) 
(30c) 
(30d) 



where the critical wave number k c , and the periodic 
eigenfunction p\(t) are determined from the linear sta- 
bility analysis. Specifically, we take (fc c ,/ c ) to be the 
minimum of the neutral stability curves and use i|10fl . 
evaluated at k = k c and with |^4/- C! „| = 1, to obtain pi(t). 
The phase of Ak a , n is determined by the map l|19fl . evalu- 
ated at (k, /) = (fc c , f c ), and the (subharmonic) condition 
pi(t + 2tt) = -pi(t) (equivalently, A kct1l+2 = -A k ^ n ): 



Al c ,n _ cos(27ro; fcc ) - F n c n s n+1 
A r k sin(27rcj fec ) - F n s n s n+1 



(31) 



This expression defines the critical eigenmode associ- 
ated with the instability. The subharmonic assumption 
Ak c ,n+2 = —Ak. ci n is valid when the instability is associ- 
ated with the first resonance tongue, as it is for all pa- 
rameters we have investigated. Plugging (|30c|) into ltTa|) 
yields qi(t) — -r-(p'i(t) + r yk 2 pi(t)) at leading order in e. 
Since in this section k will generally be fixed at its crit- 
ical value, we hereafter drop this subscript from many 
expressions (e.g., we write A n for ^4fe ci „) or replace it 
with a subscript 1 (such as with pi(t)). 

At second order we find the equation (cf. JHJ) 

q\ = p' 2 ' + 2 1 (2k c ) 2 p> 2 

+ h 2 (2k c f + r (2fc c ) 3 + (Go - G(t))2k c ] P2 , (32) 

which has the solution 

p 2 (t) = a ^2 e 2(- 7 fc2 +Ml)(t _ t „) + l fe | Ai |2 e -2 7 fc^- t „) 



+ B n e { ~ l{ - 2ka)2+iLJ2) ^ tn) + 



c.c, 



(33) 
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for t S (t n ,t n+1 ). Here 



-2k c Lof 



wf - 4cj 2 + 47 2 fc2 + 8i7fc 2 wi ' 



b = 



(34) 



(35) 



and loi (L02) is the frequency of a wave with wavenum- 
ber k c (2k c ) obtained from the dispersion relation 
Note that the "homogeneous" solution to lj32l) must be 
included (i.e., B n ^ 0) to ensure that j?2(i) is continuous. 
Using this continuity condition and the jump condition 
on p' 2 (obtained as in the linear stability calculation), we 
relate B n +i to B n through the map 

f ?+| ) = e- 47fc ' d "M 2 ,„ (^A + R n + bS n . (36) 

Here M 2 ^ n is the matrix Mk, n of (|14() with k — 2k c ; the 
vectors R n and S n are given in the Appendix. To obtain 
B n , we require that B n +2 = B n . (This harmonic condi- 
tion is due to the fact that q\, the driving term in (|32H . 
is 27r-periodic.) 

At third order we obtain the solvability condition 
(see |12) 



dZ 

t— = Lf 2 Z + {C res + C non )\Z\ 2 Z, 
where r, L, C re s and C non are given by 



1 

2^ 



47T" 



(pi +jklp 1 )p 1 dt 



(37) 



(38a) 



2tt 



[AoAoie^-l-iu!) 



k ,>47r 

4tt j . 



+ Aili(e iwidl - 1 - iw x )) + c.c, 

Gimp (^) 



/ 



-PiPidt 



(38b) 



(4,4, + 44 - 44 - 44) + c.c. 

.2 /^Tr" 



2n 

"c_ 

2- 

fc 3 



[(91P2)' + jklqip 2 ]pidt, 



(38c) 



0- 



C no „ = ^ / [-(Pi9i)' - 7*£Pi9i + k c qtPi (38d) 
3 

+ ^Top^pidt. 

The full expressions for C res and C rafm are provided in 
the Appendix. In equations we use pi to indicate 
the solution of the adjoint linear problem, which has the 
same form as ||SJ| but with 7 — > —7. Hence 



Pi 



(t) = ^etTk'-Ha.OCt-tn) + c . c . ; t e (t n ,t n+1 ) , (39) 



with the complex conjugate of A n denoted by A n . The 
map relating 4i+i to A n is similarly obtained from i|13|) 



by replacing 7 with —7 and setting k = k c . Note that in 
(|37|l we separate the cubic coefficient C into two distinct 
contributions, C res (resonant) and C non (nonresonant). 
The C res term can be traced to quadratic nonlinearities 
in Q (i.e., terms in which the k c and 2k c modes interact), 
while C non derives from cubic nonlinearities and involves 
only the critical k c mode. 



B. Weakly nonlinear results 

We now examine in greater detail the cubic coefficient 
C = C res + C noni which determines the nonlinear satu- 
ration of the instability to standing waves. Figs. BHa,b) 
show C as a function of the capillary parameter To for 
various choices of A. We find, for most parameters, that 
C(Fq) < 0, ensuring that the bifurcation to standing 
waves is supercritical. (We find subcritical bifurcations 
only in the capillary wave regime, To ~ 0.25, and with 
extreme asymmetry, e.g. A w 0.99.) 

A dominant feature of the curves shown in Figs.|3a,b) 
is the dip in C that reaches a minimum value at To = 
r re s ~ 0.09. This resonance feature is apparent both 
with impulsive forcing (non-solid lines in Figs. Eta,b)) 
and with sinusoidal forcing (solid lines). For this figure 
we have used a normalization convention for p\ (t) in the 
sinusoidal case which agrees with our normalization con- 
vention in the impulsive case (|4 c ,n| = 1) in the limit 
7^0. Note that the sinusoidal and impulsive results 
agree quantitatively for 7 = 0.01 and A = 0. For larger 
damping (7 = 0.05), the sinusoidal curve is shifted rela- 
tive to the impulsive ones indicating a difference between 
their nonresonant contributions C non to C . 

In the sinusoidal case, the dip in C is easily attributed 
to a 1:2 spatio-temporal resonance, i.e., at r res we have 
2uji = 0J2 (see 0, for example). The strength of this 
feature, which scales with I/7, is due to the strong cou- 
pling between the first subharmonic and the first har- 
monic modes, which become spatially commensurate at 
r res . Since harmonic resonance tongues are entirely ab- 
sent from the neutral stability picture for equally-spaced 
impulses, it is perhaps surprising that the dip appears in 
this case as well. However, we can understand its origin 
in the following way (see [32| for an analogous argument 
in the case of two-frequency forced Faraday waves). We 
express the height of spatially-periodic surface waves at 
integer multiples of the forcing period in the form 



H 1- c.c, 

where t m = 2ixm, m G Z. Here fi is the (real) Floquet 
multiplier of the k c mode while A and A are the Flo- 
quet multipliers of the 2k c mode(s), assumed to be com- 
plex. From the linear stability analysis, we know that 
the critical standing waves are subharmonic, i.e., at on- 
set fi = — 1. Considering the symmetries of the problem, 
including spatial translation and reflection, we arrive at 
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FIG. 9: Cubic coefficient in §ffl\ as a function of To with (a) 
7 = 0.01 and (b) 7 = 0.05. The dashed line indicates the 
result for impulsive forcing with A = 0.5; the dot-dashed line 
is for A = 0; the dotted is for A = —0.5; the solid line is for 
sinusoidal forcing, (c) 1 — A as a function of To where A and 
A are the Floquet multipliers I2H of the 2k c mode for A = 
and 7 = 0.01 (solid line) or 7 = 0.05 (dashed line). 



a stroboscopic map from the m th to (m + l) s * period: 



z m +i = ~z m + {qiWm + q.\V m )z m + a\\z m \ z m + 
2 1 



w m+1 = Xw 
v m+1 = Xv r 



~ 12 Z r: 
92 4 



(41) 



where we have kept the terms that contribute at leading 
order to the cubic coefficient in IpTfjl . Because w m and 
v m are damped modes (|A| < 1), we may perform a cen- 
ter manifold reduction to eliminate them. This gives, at 



leading order, 



92 



(1-A) 
resulting in the map 



12_ 2 

— . 



(1-A) 



-z m + \ 01 



<n<i2 



1l<l2 



(1-A) (i-A) 



(42) 



(43) 



If A and A are close enough to one, then the slaved modes 
can make a large contribution to the cubic coefficient C . 
Using the Floquet analysis of Section III Bl we confirm 
that for r = r res the Floquet multipliers associated 
with the 2k c mode come closest to one, i.e., at that point 
the Floquet multipliers cross the positive real axis (see 
Fig. Etc)). This resonant contribution to the cubic coeffi- 
cient C becomes increasingly pronounced as the damping 
decreases (compare the vertical scales in Figs. Eta) and 
(b))- 

A comparison of the (non-solid) curves in Figs. EIa,b) 
reveals that the spacing of the impulses can have a sig- 
nificant effect on the magnitude of the resonance feature. 
Some explanation for this behavior is provided by recent 
results [3, H| on multifrequency forced Faraday waves in 
which the form of the cubic amplitude equations is ob- 
tained from the spatial symmetries and the (weakly bro- 
ken) symmetries of time translation, time reversal, and 
Hamiltonian structure. Such a procedure yields impor- 
tant qualitative features of the behavior of the critical 
modes and other damped modes driven through non- 
linear interactions, including a prediction about which 
particular damped modes are most important and which 
forcing frequencies (and corresponding temporal phases) 
are most effective in enhancing or otherwise controlling 
their effect. Here we are interested in the Fourier se- 
ries expansion for an impulsive forcing function with two 
impulses per period: 



Gimp (i) 



f. 



imp V^/i 
271" ^ [ 



e -»i((l+A)7rh e »it 



+ c.c. 



(44) 



I n 0- Q the contribution of the damped 2k c mode to 
the 1:2 spatio-temporal resonance is found to be most 
affected by the forcing frequency 2lu (where ui is the pri- 
mary driving frequency) and its phase. This result sug- 
gests that we focus on the drastically truncated two-term 
Fourier series 

jimp IY-, , _ ijrA 



G 2 (t) 



2tt 



[(1 



)e lt + (1 - e - 2lwA )e 2lt + c.c 



(45) 

At leading order in the damping 7, the onset of standing 
waves occurs when the magnitude of the first term in 145|) 
becomes equal to 7 (see, for example, 0). Using 



imp 



1 



^ — in A I 



(46) 



and simplifying the result through an appropriate time 
translation, we may rewrite l)45|) at onset as 



G 2 {t) = 7 e lt + F 2 e 2lt + c.c. 



(47) 
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where 



(a) 



F 2 = 2i7sin (A- t 



(48) 



From Table 1 in Q we find the predicted magnitude, 
at leading order in 7, of the resonant dip in the cubic 
coefficient at r rP .,, 



C 2f . 



-ai- 



\L 3 \ +Mi|/2a|sin$ 



(49) 



Here |/ 2a | = |F 2 | = 2 7 sm(A§), $ = arg(F 2 ) = 
sign(A)?-, L3 is the linear damping coefficient of the 2k c 
mode (L3 = (2fc c ) 2 7 w 47), and /Zj is the coefhcient 
of the linear parametric driving term for the 2k c mode 
(fJ-i = 2 ~ !)■ Thus Cy/ a , as a function of the asymme- 
try parameter A, becomes 



a 2 / s (A)«- 



27 2 



1 



■(¥) 



(50) 



The constant ot\ is positive and must be determined by a 
nonlinear calculation; it depends on To, but not, at lead- 
ing order, on 7. Here we choose ol\ such that the values 
of the multifrequency prediction Cy/ S (A) and our direct 
calculation of C res (A) at Tq — T res from l|A3|l agree for 
small damping (7 = 0.001) at A = 0. (We then use 
the same estimate of u\ throughout this comparison.) 
We compare C res (A) with the multifrequency prediction 
in Fig. 1101 and find good agreement, especially for small 
7 where the multifrequency results are expected to be 
valid. This agreement breaks down as |A| approaches 
unity (even more so for larger 7), an effect that is not 
overly surprising given that |A| — > 1 is an unphysical 
limit requiring an infinite value of fi mp to produce stand- 
ing waves. Note that C%£ s , given by l|50|) . scales with 



7" 



which explains the difference in scales evident in 



Fig.EKb). 




IV. DISCUSSION AND CONCLUSIONS 

In this paper we examined the problem of Faraday 
waves parametrically excited by a periodic sequence of 
impulses (delta functions). We showed how to extend 
the linear stability analysis presented in 0] to include 
more general impulsive forcing functions, and investi- 
gated the weakly nonlinear regime for one-dimensional 
surface waves in the case of two impulses per period, com- 
paring both linear and nonlinear results with the more 
established cases of sinusoidal and multifrequency forc- 
ing. Our analytical and numerical work was conducted 
within the framework of the Zhang- Vinals Faraday wave 
model, appropriate for describing small amplitude sur- 
face waves on deep, weakly- viscous fluid layers. 

One theoretical advantage of using an idealized im- 
pulsive forcing function is that it allows for an exact, al- 
beit generally implicit, expression for the neutral stability 
curves associated with the primary instability. Moreover, 



FIG. 10: Resonant contribution C re s to the cubic coefficient 
of Eq. 1)3 7fl as a function of the asymmetry parameter A over 
the range [—0.95, 0.95], calculated for the Zhang- Vinals equa- 
tions using l|A3ft (solid line) and using multifrequency forcing 
results 150H with a\ — 0.235 (dashed line). Plots are for (a) 
7 = 0.001, (b) 7 = 0.01, and (c) 7 = 0.05. 



since the stroboscopic map characterizing the linear prob- 
lem can be explicitly constructed, Floquet multipliers can 
be readily determined as a function of forcing and fluid 
parameters. This stands in sharp contrast to the sinu- 
soidal and multifrequency cases where even the neutral 
curves must be determined numerically or through an 
asymptotic expansion that assumes weak damping and 
forcing. A further consequence of exactly solving the lin- 
ear problem is that we arc then able to derive explicit, 
analytic expressions for the amplitude of weakly nonlin- 
ear, one-dimensional surface waves as a function of forc- 
ing and fluid parameters. 
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As noted by 0] , in the simplest case of TV = 2 impulses 
per period, harmonic instabilities from the flat fluid state 
cannot be excited with equally-spaced impulses. Despite 
this fact, we found, through a weakly nonlinear analysis, 
that the 1:2 spatio-temporal resonance remains the dom- 
inant feature in typical parameter regimes. Symmetry 
arguments of the type used in [3^ showed that, as in the 
case of sinusoidal forcing, harmonic modes can be driven 
nonlinearly by the critical modes — hence the resonance 
effect persists. 

By varying the spacing between the up and down im- 
pulses making up the 27r-periodic forcing function for 
N = 2, we found that the magnitude of the 1:2 resonance 
effect depends dramatically, and monotonically, on the 
corresponding asymmetry parameter A. Appealing to re- 
cent results from multifrequency forcing theory [4|, valid 
in the limit of weak damping, we obtained a prediction of 
this dependence based on a truncated two-term Fourier 
series approximation. This prediction agreed quite well 
with the results for impulsive forcing despite the severity 
of the approximation involved. 

We envision several ways in which this work could be 
extended. For example, the methods used in this pa- 
per could be easily applied to the case of piecewise con- 
stant forcing, a type of driving function that has been ex- 
plored in other systems (see [3(| |3^ , for instance) where 
it also led to analytic results. It would further be of 
interest to extend our weakly nonlinear analysis to the 

I 



and 



The expressions for a and b in these vectors are given 
by (|34|l and (|35|l . respectively, and u>i (0J2) satisfies the 
dispersion relation Q with wave number k c (2fc c ). 



case of two-dimensional patterns. In such a context, ad- 
ditional questions of pattern selection, as well as more 
specific comparisons with multifrequency forcing, could 
be explored. One could try, for example, to construct 
impulsive forcing functions that mimic specific effects 
seen with multifrequency forcing. In this way the an- 
alytic results available for impulsive forcing would com- 
plement the numerical I23l I26L l33l and experimental re- 
sults ElEQHHTIBUliil reported with multi- 
frequency forcing. 
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APPENDIX A: ANALYTIC EXPRESSIONS 

The vectors given in l|36|) are 



(Al) 



(A2) 



I 

The full expressions for the components (C res and 
Cnon) of the cubic coefficient in the solvability condition 
iJSZJ are 



Rn = 



Re[a(A 2 n e ( --^+^ 2d - - A 2 n+1 )} 



5,i 



i(K| 2 e- 2 ^-K +1 | 2 ) 
-27*X((_i)«2k c / c - 7 fc c 2 ) + K +1 | 2 7 fc c 2 ) 
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C res — 



Cn 



E AjBA&^ief** - 1) - A^^p^ie^ - 1) (A3) 
Pi Pi 



i=o,i 



+ (Wi^ - A jBj A/-^-) (e** - 1) 

+ (aA^^l + {h -a)\A^A A 3 ^l) (e** - 1) 



+aA^A j ^l(e^ - 1) + (6 - a)|^f ^^^(e" 2 ^ _ 1} + c . c ., 



E + 3r ° fc %;; ^V *" - 1) (A4) 



j=0,l 



where 

/3± = -4 7 fc2 + i(±2wi + w 2 ), 02 = -47^c + ^2, 03 = -2lk 2 c + 2iwi, /3 4 = -2 7 fc2 + 4iwi. (A5) 
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